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O^' We propose an indirect method for observing radiation from an incomplete soliton with a suffi- 

"t^ ' ciently large amplitude. We show that the radiation causes a notched structure on the envelope of 

the wave packet in the momentum space. The origin of this structure is the interference between the 
i-A ' main body of oscillating solitons and the small radiation in the momentum space. We numerically 

G ■ integrate the nonlinear Schrodinger equation and perform Fourier transformation to confirm that the 

Q_) ' predicted structure really appears. We also show a simple model which reproduces the qualitative 



result. The experimental detection of the notched structure with the Bose-Einstein condensation of 
neutral atoms is discussed and suitable parameters for this detection experiment are shown. 
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1. Introduction 

^ ! The one-dimensional nonlinear Schrodinger equation (NLSE), 

H ■ 

i«Ar - -2'A« - l«AlV, (1) 

is a well-known soliton equation which appears in various fields of physics^ and has many applica- 
tions, such as in optical fiber communication.^'^ Another important example of physical systems for 
the NLSE is seen in the Bose-Einstein condensation (BEC) of neutral atoms,"*' ^ where the macroscopic 
wave function of condensate atoms appears as the order parameter accompanied by the spontaneous 
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breakdown of the L'^(l) gauge symmetry.^ In this case, the NLSE is regarded as the mean-field ap- 
proximation of the Heisenberg equation for field operators and describes the time evolution of this 
macroscopic wave function with good accuracy.'*'^ Up to this day, both bright and dark soliton-like 
objects have already been created with BEC^'^ (though the dark solitons are outside the scope of this 
study.). 

On the other hand, it is clear that we cannot realize a rigorous soliton initial condition in any 
real experiment. However carefully prepared, any realized physical object must be accompanied by 
deviation from ideal solitons. Even in such cases, it is widely known that the NLSE has analytical 
integrability; exact solutions under an arbitrary initial condition are made available by the inverse 
scattering transformation (1ST) method.^' ^^ In particular, to analyze soliton phenomena, we can exploit 
Hirota's direct method. ^ ' The 1ST method gives us not only the solutions of the NLSE but also a clear 
outlook on their classifications. It is based on an auxiliary linear eigenvalue problem. 



»Pij_ -i^ ir ppi 



(2) 



where ^ is the eigenvalue which is independent of time, ip is the solution of the NLSE, and the asterisk 
means its complex conjugate. This formulation takes the form of a potential scattering problem for the 
auxiliary fields ^i and ^2. where ip works as a potential. The eigenvalue spectrum consists of discrete 
and continuous parts, where the former generates soliton solutions and the latter corresponds to small 
ripples called "radiation". For a sech-type initial condition. 



ipix, 0) = Msech(x), 



(3) 



the eigenvalue consists of a discrete part only, provided that M equals a positive integer. In this case, 
the whole initial value problem is solved analytically and results in an M-soliton bound state. ^^ 
The time evolution of the auxiliary field is defined by another linear equation. 









l^2J 



(4) 



i^^ - iii«Ai' 5<a; - i^r 

-i^,-i^«A -i^" + ^i|«Al^ 
However, the time evolution of the wave packet from initial conditions except eq. (3) is relatively 
unclear since the analytic treatment of the continuous spectrum is difficult. All we know is the asymp- 
totic behavior. According to the analytic results obtained by Satsuma and Yajima,^^ a solitary pulse 
is very robust and, unless the spectrum of the eigenvalue problem (2) contains a discrete eigenvalue, 
the wave packet finally splits into the soliton and radiation parts. The norm of the soliton part remains 
0(1). When the amplitude of the wave packet is sufficientry large, the bound soliton can be formed 
asymptotically and has already been observed in real experiments with optical fibers. ^^ In contrast, 
the amplitude of the radiation part decays as 0{r^) and expands rapidly from the main body of the 



2/16 



J. Phys. Soc. Jpn. Full Paper 

soliton part. 

In this way, the amplitude of radiation is generally very small, which makes it difficult to observe 
radiation in real experiments. In fact, whatever the research field is, no clear observation of radiation 
has been reported. The study of radiation has been difficult both theoretically and experimentally, and 
the problems on the dynamics of the NLSE, starting from general initial conditions, have not been 
solved. 

In this paper, we propose an idea toward the experimental observation of radiation through pat- 
tern formation in momentum space. We exploit the function of order formation originating from the 
nonlinearity of the NLSE. 

This paper is organized as follows. In the next section, our theoretical idea is explained and the 
results of numerical simulations are shown. In §3, the feasibility of a real experiment is discussed. 
Section 4 is devoted to discussion and summary. 

2. Theory and Numerical Results 




Fig. 1. Profiles of the initial wave packets. The dashed curve is for intial condition (5) and the solid curve is 
for initial condition (7). 



In this section, we briefly summarize the mathematical descriptions of the system to be considered. 
Throughout this paper, we consider the NLSE (1) under initial conditions, the eigenvalue spectrum of 
which contains two or more discrete eigenvalues. By virtue of the Galilei transformation, we can 
restrict ourselves to the case where the center of the wave packet remains at x = 0. Such a situation 
is achieved as long as the initial conditions are real valued and symmetric functions of x. Under such 
conditions, the relative velocities of the contained solitons vanish and these solitons form bound states. 
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Fig. 2. Momentum space profile of the wave packet \\p\{k, f)P starting from initial condition (5) at t - 14. 




Fig. 3. Momentum space profile of the wave packet \ij/2{k, 01^ starting from initial condition (7) at f = 14. The 
dashed curve is for the initial profile. 

We consider two similar initial conditions which meet the requirement above. The first initial 
condition is the sech-type (3) one with M = 2, 



iAi(x,0) ^2sech(x). (5) 

Starting from ij/i{x, 0) results in a pure two-soliton bound state, and the solution of which reads, 

it cosh(3x) + 3exp(-4i?)cosh(x) 
(Ai(x,0 ^ 4exp(--) , ,, , , ,__,,..^ ^ . o__,...x - (6) 



2 cosh(4x) + 4cosh(2x) + 3cos(40 
The envelope of this solution pulsates with the frequency | and no radiation is emitted. 
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Fig. 4. Momentumspaceprofileof the wave packet It/rsoi+tAradl, where /(A:) is taken to be e"2 . This snapshot 
was taken at f = 14. 




Fig. 5. Profile of the ground state wave packet in the external potential Vf^y^xix) - \x^. The real line is for the 
profile of the wave packet and the dashed line is for the external potential yext(''^)- 



The second initial conditon is of the Gaussian type with the same amplitude and a similar width, 

1 



tp2{x,Q) = lex^{--x ). 



(7) 



This form of initial condition not only deviates from the complete soliton but also becomes one of the 
accessible forms, especially in BEC experiments, and its time evolution was investigated in detail in 
our previous work.^^ By the Feshbach resonance technique, ^^ we can erase the effective interatomic 
interaction, and the atoms trapped in a quadratic potential trap take the ground state to form the Gaus- 
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Fig. 6. Momentum space profile of the wave packet starting from the ground state one in the external potential 
Vextix) = \x^. This snapshot was taken at f = 14. The dashed curve is for the initial profile. 




Fig. 7. Profile of the ground-state wave packet in the external potential Vextix) = ^x^- The real line is for the 
profile of the wave packet for g'j^ = 0.7, and the dashed line is for the external potential yext(-*^)- 



sian profile. We show |(Ai(x, 0)P and \tfr2{x,0)\^ in Fig. 1. They show similar curves; in fact, both lead 
to pulsars. However, the wave packet starting from initial condition (7) does not have a simple ana- 
lytical expression and keeps emitting a very small amount of radiation which has not been observed 
experimentally. 

We propose a method of distinguishing these two initial conditions. 

For the above purpose, we observe that there appears a distinguishable feature on the envelope of 
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Fig. 8 . Momentum space profile of the wave packet starting from the ground-state one in the external potential 
Vextix) - \x^. This snapshot was taken at f = 14 with gj^ - 0.7. The dashed curve is for the initial profile. 
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Fig. 9. Profile of the ground-state wave packet in the external potential yext(-*^) - \x^- The solid line is for the 
profile of the wave packet for g'^^ = 5, and the dashed line is for the external potential VextC-^)- 



a wave packet which emits radiation. That is, a "notched structure" is formed on the envelope in the 
momentum space. The creation of the notched structure is explained as follows. 

When we consider an initial condition with a sufficiently large amplitude that deviates from that 
of the complete soliton, the total wave function can be expressed as the sum of a soliton part i/'soi and 
a radiation part ij/radL, i-c, lAsoi + (Arad- We write these two constituent wave functions as (Asoi = g(x)e'^'^^^ 
and i/'rad - /j(x)e"'^^\ where g{x) and h{x) are absolute values of tAsoi and t/^rad, respectively. They have 
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Fig. 10. Momentum space profile of the wave packet starting from the ground-state one in the external po- 
tential Vf^xiix) - \x^. This snapshot was taken at f = 14 with g'^j^ - 5. The dashed curve is for the initial 
profile. 

a definite pliase 6{x) in common; tlius, they never interfere with each other in the real space. However, 
when we perform the Fourier transformation on these wave functions, the difference between g{x) and 
h{x) induces two different phases for i^goi and i^rad» where the tilde means a Fourier component. Now, 
they can interfere with each other in the momentum space. This interference pattern in the momentum 
space is nothing else but our notched structure mentioned above and this is an outline of our strategy. 
Below, we show a brief calculation of the expression of the interference pattern on the bases of a rough 
but concise model. 

Since the pulsating soliton part localizes around the origin and oscillates with monochromatic 
frequency oj = |, it can be approximated as 

lAsoiC-^, ^ e^5^(3 -I- cos(40). (8) 

We note that this approximation holds only around the origin; in fact, when ifj{x, t) shrinks, the rest 
part of the main body of solitons extends toward the left and right to conserve the whole norm. How- 
ever, this rather rough approximation does not spoil the qualitative argument because the oscillating 
amplitude of the side lobe is much smaller than that of the main body of solitons near the origin. 

On the other hand, we can neglect the nonlinear interaction term of the radiation part because 
its amplitude is very small. Therefore, the radiation part is well described by the linear Schrodinger 
equation. Such a wave motion with zero group velocity is expressed as 

"AradCx, - J fik)c-i'\'''''dk, (9) 
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where f{k) is a real and symmetric function of ^.This expression for the radiation part is no more 
than an approximation. In fact, there exists an interaction between the radiation part and the soliton 
part, particularly near the origin where the amplitudes of these parts are large. In such an area, the 
approximation is poor. However, such an area is relatively very narrow and limited to around the 
origin. To observe the notched structure, we perform the Fourier transformation on the wave function 
in the real space. By virtue of Fourier transformation, which is an integral over an infinitely long 
interval, the local error is smeared or averaged out in the momentum space and brings about no crucial 
alternation but a small offset to our original idea. In almost all othe areas but the origin, we have 
confirmed numerically that our approximation is valid. 

We shall consider the interfered profile of these soliton and radiation parts in momentum space, 

l^sol + ^radP. (10) 

For the validity of superposing these approximated functions, see the appendix. Obviously, we find 

^soi(^, = e"5*^'(3 + cos(40) and iArad(^, = /We"^*^''. Therefore, 

l«Asoi + «Aradl' - (3 + cos(40)2e-'=' + p{k) + 2(3 + cos{4t))mc--2''\os{jk^t). (1 1) 

This interference pattern is observed as the notched structure on the surface profile of the wave packet 
in momentum space. The larger k makes the oscillation period shorter and the structure is sent out of 
the center of the profile. This spatial and temporal order is quite outstanding and it makes it possible 
for us to observe radiation indirectly. 

We show the results of numerical simulation based on the model explained above. 

In Figs. 2 and 3, the snapshots of \tl/\{k,t)\^ and 10-2(^,01^ are depicted at f = 14. Although the 
surface of \i^i{k, t)\^ is smooth, the notched structure is fully developed on the surface of li/'iC^, OP- For 
reference. Fig. 4 shows |iAsoi + i/'radl att - 14, where f{k) is taken to be e^2 . It is seen that our model 
qualitatively well describes the result in Fig. 3 even though the norm becomes much larger because 
our simple model does not preserve the norm. 

Next, we consider a more general case. In real experiments, one usually trap atoms with an external 
potential Ve^i{x), which can be approximated by a quadratic function. During the trapping time, the 
macroscopic wave function of atoms goes to the ground-state. The ground-state wave function in the 
potential corresponds to the initial condition. After this process in the trapping potential, we switch 
off the trap and release the atoms to observe the interference pattern. Figure 5 shows the profile of 
the ground-state wave packet in the external potential KxtC-^) - 5-^^- The real line is for the profile 
of the wave packet and the dashed line is for the external potential. Figure 6 shows the corresponding 
momentum space profile of the wave packet starting from the ground-state one in the external potential 
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Vext(-'c) - jx^- This snapshot was taken at t = 14, which is sufficientry long for the wave packets to 
develop the interference pattern. Figure 6 clearly shows the emergence of the notched structure from 
a more general initial condition than that in a realistic experimental setup. 

3. Experimental Feasibility 



In this section, we make some remarks on the possibility of real experiments with the BEC of 
neutral atoms. One of the striking features of the BEC experiments is the visibility of the momentum 
space. One can extract the momentum distribution of cold atoms by the time-of-flight (TOE) method. ^^ 
In the TOE method, while atoms with large momenta can travel far away, atoms with small momenta 
can travel only small distances from the origin. Therefore, the momentum distribution is detected 
directly and we can observe the small radiation ripple through the formation of the interference pattern 
with the aid of the TOE method. Moreover, the control of external environments is relatively easy in the 
BEC systems where we can confine condensate particles along a quasi-rectilinear line by tightening the 
laser beam trap. Therefore, the confirmation of the emitted radiation in BEC experiments is promising. 
As general requirements for this experiment, we note that a sufficiently large amplitude is required to 
form the 2-soliton bound state and the strength of the external potential is important. If one uses a very 
loose potential, the ground-state wave packet might have an insufficient amplitude and might fail to 
form the 2-soliton bound state. We recommend that the coefficient in front of x^ be chosen as larger 
ffian J. 

In addition, we must take special care for some subtleties arising from one-dimensional geome- 
try and attractive interaction. In a recent study, Tovbis and Hoefer^^ also treated the self-interaction 
and pattern formation of solitons in a similar environment and discussed the possible experimental 
condition. Although their argument is for the real space, which is different from our argument in the 
momentum space, we should also make considerations on the experimental feasibility. Eor ffiat pur- 
pose, we rewrite the NLSE (1) back in the original form with dimensional constants as 

in-^ = -—-^-gwm^^, (12) 

ot 2m ox^ 

where h is the Dirac constant, m is the mass of the Bose-condensed atoms, and gio is the coupling 
constant. In the typical three-dimensional geometry, ^30 is written as g^o = " , where a^o < is 
ffie s-wave scattering length of the atoms with an attractive interaction. The s-wave scattering length is 
a quantity peculiar to a three-dimensional geometry and it is not a trivial problem to seek for the one- 
dimensional correspondent. However, in real experiments, we can stay in the NLSE regime provided 
that the radial confinement diameter aj_ is experimentally suitable level. ^^ 

The next issue is related to the realization of eq. (1). In a typical experiment, we initially trap atoms 
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by applying a harmonic external field and suddenly switch off the trap and release atoms to recover 
the fiat geometry along the axial direction at t = 0. Strictly speaking, when we use an optical trap, a 
relatively loose Gaussian potential remains along the axial direction even after the release. However, 
its curvature can be neglected for our purpose. The harmonic external field defines the typical length 
scale of the system, Uho = J-^^ where co is the trap frequency. We can reduce eq. (12) to eq. (1) by 
replacing x with J—x, t with -, and >fj with J^i^if', where A'^ is the total number of particles, and 
by letting Jl° be unity. It follows that 

A^^i:^. (13) 

gw 

More generally, if we let ^J.° be .gj^,, a dimensionless coupling constant, eq. (1) is generalized to be 

1 

2 

and we obtain 



iil^t = --^xx - g'w\^\^^, (14) 



A^=— -^. (15) 

g\D 

In real experiments, it would be difficult to adjust g'^^ exactly to unity. For such cases, we give here 
the tolerance for g\j~^ to detect the desired fringe. Figure 7 shows the profile of the ground-state wave 
packet in the external potential Vf^xt{x) = ^x^ with g^^ kept at 0.7, and Fig. 8 shows the corresponding 
momentum space profile of the wave packet starting from the ground-state one in the external potential. 
This is a snapshot taken at t = 14. It is found that the contrast of the fringe is relatively low because 
a smaller coupling constant g'^^ makes the role of nonlinear oscillation less important. Therefore, we 
can roughly estimate that the lower bound of the suitable g'^j^ is about unity. 

Next is the upper bound. Figure 9 shows the profile of the ground-state wave packet in the external 
potential yext(-'c) - \x^ with gj^ kept at 5, and Fig. 10 shows the corresponding momentum space 
profile of the wave packet starting from the ground-state one in the external potential. This also shows 
a snapshot taken at t - 14. In this case, a strong self-focusing interaction makes the initial profile 
narrower and this leads to a wide-spread profile in the momentum space. In this case, the image 
intensity in the momentum space would be lower. We should keep image intensity above the detectable 
level, and it defines the upper bound of g'^j^. As a whole, the suitable range for the fringe-detection 
experiment is roughly estimated to be about 1 < g'^j^ < 5. 

Finally, let us consider preparing initial condition (7). As previously mentioned, we can change 
g\D freely by the Feshbach resonance technique to gain a sufficient number of atoms, as seen from 
eqs. (13) and (15). However, the order of A'^ is upper-bounded, because this quantity must satisfy the 
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inequality 

N < 0.67-^, (16) 

to avoid the self-collapse of the condensate owing to the attractive interaction. ^^ In real experiments, 
Strecker et al. used the values a^ k 1.5 /um and a^D k -0.16 nm, letting the order of N be several 
thousand. Therefore, the experimental feasibility of our proposal seems to be promising with the aid 
of the Feshbach resonance technique, provided that the number of particles is kept moderate. 

4. Discussion and Summary 



As evidence that supports our simplified model, we can consider that our model explains why "the 
higher the wave number, the narrower the pitches of the interference pattern," as shown in Fig. 10. 
This is due to the term including cos^k^t in eq. (11), which validates the low-amplitude-defusing 
wave approximation of eq. (9). 

Next, the oscillations in the momentum space itself are not rare phenomena. As an example, 
we can immediately give the rectangular function which leads to the sine function after the Fourier 
transformation. However, such functions tend to be angular, in general, and to exhibit oscillations in 
the momentum space. On the contrary, smooth functions, such as the Gaussian and sech-type ones, do 
not exhibit oscillations even after the Fourier transformation. Our wave packets which smoothly and 
monotonically decrease from the origin are among them and are not angular. Therefore, "interference" 
by our model is the most promising mechanism to explain the oscillation in the momentum space. 

We have proposed an indirect method for observing radiation from an incomplete soliton. Accord- 
ing to this method, radiation forms the notched structure on the surface profile of the wave packet in 
the momentum space with the aid of the TOF method. The origin of this structure is the result of the 
interference between the main body of the oscillating soliton and the small radiation in the momentum 
space. We numerically integrated the NLSE and Fourier-transformed it to confirm that the predicted 
structure really appears. Our concept was successfully proved and our simple model reproduced the 
qualitative result from the exact numerical simulation. Finally, we discussed the possible conditions 
required in real experiments with BEC. 
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Appendix 

In this work, we consider initial conditions that sUghtly deviate from that of the exact 2-soUton 
bound state (5). Their eigenvalue spectra should include 2 discrete eigenvalues. The typical example 
is the Gaussian wave packet (7). In fact, initial conditions (5) and (7) appear very similar, as shown 
in Fig. 1. The stationary solutions formed in the harmonic potentials are also considered to meet the 
above conditions. 

Initial condition (5) evolves into eq. (6). The purpose of this appendix is to explain the validity 
of the assumption of superposition we made in eq. (10). Under this assumption, general wave packets 
that meet the above conditions can be approximated as the sum of the oscillating main part and the 
small diffusing part. The smallness of the latter enables us to express it as a general solution of the 
linear Schrodinger equation. 

Consider the linear potential scattering problem (2) and give the corresponding lost functions as 



<p{x;^) = 



( \ 
1 



exp(-/^x), X -^ -oo. 



«A(x;^) = 







exp(/^x),x -^ oo. 



4f{x;^) = 



1 







exp(-/^x), X ^> oo., 



where £, is the eigenvalue of the linear scattering problem (2) and consists of discrete and continuous 
parts. We write <p{x\ ^) as a linear combination of ^{x; ^) and i/r(x; ^), 

0(x;^) = a(^)«A(x;^) + ^(^)^(x;^), 

where the transmission and reflection coefficients a(^) and h{^) satisfy 

\a{^f + mf = 1- 

From the general arguments of the potential scattering theory, ^ can be analytically continued to the 
upper complex plane, and bound states correspond to the discrete zero points of a{^), 

a(4) - 0. 



According to Zakharov and Shabat, the solution u{x, t) of the original NLSE (1) can be written as 






1 r°° 

exp(-/^„x + /^20<A:2 + - ^\i^)d^- 



Here, the zeros of a{^) are assumed to be simple poles and a'(^„) - ^lf=f„ is not equal to zero. In the 
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above formula, the discrete part including t/^„2 corresponds to the soli ton part {n = 1, 2 for the 2-soliton 
bound state) and the integral part including O2 to the small diffusing wave. This is clearly a superposed 
form of the soliton-like and small-diffusing-wave parts. Neither can be expressed analytically. The 
integrand <I>2(^) is the quantity that is determined by the simultaneous equations below. 



<Di-c(x,^) 



1 

2(1 + T) 



XL 



^ = -cix,oY,TrF<2' 



^-Cn 



c*{x,0 



2(1 - T) 



Oi+0*-c*(x,^) 1 + 2] 



A 



^-^n 



^n\ 



^m ^ bin ~ bn ^^' J-oo s ~ btn 

- > ^— (Am + -^ - 1 + — —^, 

„ bm bn ^m •^^' ^-oo S bm 



where we introduced the abbreviations 
and 



^^(4) ,.^ .^2,x 



a(^) 



exp(2/^x - 2/^ f)- 



The Hilbert transformation operator T acts as 



rd) 



1 r~ (D(^' 



J^. 



Since the pure soliton-like eq. (5) acts as the reflectionless potential in eq. (2), h{^) should be zero. This 
makes c(x,^) = 0, and we can easily see Oi = <1)2 - from the above simultaneous equations. Then, 
the small-diffusing-wave part completely disappears and only the discrete soliton part corresponding 
to eq. (6) remains. This is a pulsar localized around the origin and oscillates with period |. 

In this paper, we consider the initial conditions that slightly deviate from that of the exact 2-soliton 
bound state (5). From the principle of continuity, they are considered to evolve very similarly to pulsars 
(exact 2-soliton bound states). Because the variation in magnitude at the origin is from 2 to 4, we have 
modeled this pulsar as eq. (8). This function is localized at the origin and oscillates there with period 
|. It satisfies the condition to be a model function of the exact 2-soliton bound state (5). By numerical 
integration, we have already confirmed that the wave packets starting from the initial conditions that 
slightly deviate from that of the exact 2-soliton bound state well retain the properties of a solitary 
wave.'^ As a soliton-like part, we can use eq. (6) instead of eq. (8). However, this function cannot be 
Fourier-transformed analytically and can be unsuitable as a brief model. 
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The amplitude due to the integral part is considered to be small because the integrand $2 is small 
as long as b{^) and c(^) are small. Therefore, this part can be well approximated by the general solution 
of the linear Schrodinger equation (9). 
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